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Abstract 

The general 2-dimensional fermion system with repulsive interactions (typified by 
the Hubbard Model) is bosonized, taking into account the finite on-shell forward scat- 
tering phase shift derived in earlier papers. By taking this phase shift into account in 
the bosonic commutation relations a consistent picture emerges showing the charge- 
spin separation and anomalous exponents of the Luttinger liquid. 
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The proper description of the effect of finite forward-scattering on-shell phase shifts on 
Fermi systems for D > 1 |]T| has been the subject of a number of papers [2-4]. The existence of 
such scattering, leading to on-shell singularities in the T-matrix, was confirmed by Fukuyama 
and Narikiyo ||, di Castro, Castellani and Metzner || and others. The discovery which 
these papers confirmed followed from the fact that when two particles are embedded in 
their respective Fermi seas, effectively all soft recoils are forbidden to them by the exclusion 
principle. Under these conditions, the logarithm of the S'-matrix for relative motion retains 
a finite eigenvalue r] Q in 2 dimensions, 

So(Q) = expi[r] (Q)} 

in the limit that the relative momentum Q = — k\ — > 0, and that both states are below the 
Fermi level. (On general grounds, it seems likely that this can happen in 3D as well, but this 
has not been proven.) It is important to be clear that S is not the conventional "T-matrix" 
defined for a hole at k scattering against another at k', whose imaginary part represents an 
incoherent decay process, vanishing at the Fermi surface; S is the on-shell scattering matrix 
of particles embedded in the Fermi sea and its phase determines the boundary condition for 
their asymptotic wave-functions at the origin of relative coordinates r — r' = 0. Addition of a 
particle modifies the wave functions of all other particles, and our endeavor is to investigate 
the consequences of this fact. 

Our initial description JF| was heuristic, merely pointing out that the effect of such a 
phase shift mimics that of a change in statistics by enforcing a partial exclusion principle 
between electrons of opposite spins. We also described the type of singular interaction which 
would give energy shifts similar to those which take place, and emphasized that these would 
"trash" Fermi liquid theory. Some of the papers have focused on this singular interaction 
(specifically Randeria, et al0 and Stamp et al. M. The treatment in terms of an interaction is 
in several respects unsatisfactory, as clarified by Baskaran[[|. But even Baskaran's discussion 
does not give us a clear insight into clean formal ways to deal with the situation. 

The problem is that the effect is best thought of as a constraint on the wave-functions, 
not as an interaction. This is most clearly seen in the Hubbard model, where the effect of 
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a strong enough repulsive potential U — > oo is to enforce a projective constraint, expressed 
as the Gutzwiller projector acting on the kinetic energy in the t — J model, for instance. 
Since the exchange term also is expressible purely in terms of projected operators, the t — J 
system is confined to the subspace defined by projected operators. 

It is worth emphasizing that renormalization group derivations of Fermi liquid theory as 
a theory of the low-energy states, such as that of Shankar, implicitly assume a free Fermion 
starting Hamiltonian. If the starting problem itself is projected onto a subspace, this property 
will remain after renormalization and F.L.T. changes into the theory we shall derive. 

In general (in 2D) the constraint appears as a phase shift, which is a boundary condition 
for the asymptotic wave function in the relative coordinates of a pair of particles. Such a 
wave function is indeterminate unless it has a boundary condition both at r — r' — > oo and at 
r — r' — > Arguments in several of the original papers show that the rest of the particles may 
be satisfactorily dealt with by taking the exclusion principle into account, and multiparticle 
encounters are not crucial. 

This local boundary condition on the asymptotic wave function at r — r' — > is a 
kinematic, rather than a dynamic, effect: there is a change in the wave functions of the 
particles, not directly in their energy. We are used to this with hard core potentials: the effect 
is best expressed as one purely on the kinetic energy, not on the potential. This kinematic 
effect dominates here because the scattering region where the potential acts is small, of order 
N' 1 compared to the asymptotic region in which the kinetic energy is modified. The way to 
make this point is that such a boundary condition can actually change the dimensionality 
of the Hilbert space of allowed wave functions. In simpler terms, such a boundary condition 
forces the wave function's nodes to shift in such a way that a particle moves into or out of 
the distant boundary, so that the same volume contains N±r]/7r particle states rather than 
N. This is what is meant by a change in the dimensionality of Hilbert space. This change of 
Hilbert space occurs in ID even as a consequence of an ordinary interaction potential (hence 
the flexibility of statistics in ID) but in all other dimensions it is distinct from the kind of 
interaction effects which can be treated perturbatively. 
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The conclusion we came to is therefore that the effect of a finite phase shift is best 
modelled as a modification of the algebra of the particles, expressed in their commutation 
relations. Projected Fermions 

( C ta)proj = Ph C ta = (1 — n i-<r) C ia (1) 

do not have the same commutation relations as ordinary Fermions, obviously, but we have 
not found the Fermion representation convenient to work with. It is much simpler to use the 
bosonized representation in terms of the Fermi surface^ fluctuations || || . 

Haldane, particularly, has emphasized that the most useful description of the dynamics 
of a Fermi system is via the operators Afcp describing the position of the Fermi surface in 
/c-space, taken to be dynamical variables, functions of a coarse-grained space and time. That 
is, he argues that Luttinger's theorem holds exactly during sufficiently long-wavelength and 
low-frequency fluctuations. (Parenthetically, even the conventional derivations of Luttinger's 
theorem [pHH depend not on the convergence of perturbation theory but merely on the as- 
sumption that excitations precisely at the Fermi surface do not decay, hence the Green's 
function is real.) We define operators 

Ak Fa (Q,r,t) 

giving the Fermi surface fluctuations of spin a at a point on the F.S. parametrized by Q, and 
at coarse-grained r, t. These are the bosonic variables: they commute for different Q and r, 
and, for non-interacting electrons, for different a . We can introduce a phase variable d a of 
the wave function at the Fermi surface, which is a function of Q,r and t, and then Ak Fa is 

Ak Fa = h n ■ V9 a (2) 

where is the local normal to the fiduciary Fermi surface. 9 and Ak F , which is equivalent to 

the particle density at Q, p (Q), are conjugate variables, and have for free fermions canonical 
1 *The bosonized version of Fermi liquid theory can be equivalently thought of as the appropriate gauge 
theory in the presence of a Fermi surface, since the bosonic variable is essentially the phase of the Fermi 
surface wave function. 
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commutation relations: 



[6,p] = m8{r-r') 8{Cl-Qf) 



(3) 



As Haldane has pointed out, this representation can be motivated by the idea of expressing 
the Fermion field in terms of two real operators p and 9 



rather than by the earlier "Tomonaga" definition of p(q) as a density of Fermions J2k c t+q c k- 
This latter representation is not possible when the Fermions are projected operators. But 
we can still speak of a Fermi surface and a Fermi surface phase for each spin which satisfies 
Luttinger's theorem, hence determines the density of particles at each point on the Fermi 
surface. In this transcription of the original idea of bosonization we follow Khveshchenko 
llj- If a Fermi surface exists this implies zero-frequency modes at each point on it, hence 
separate, independent conservation of particle and spin currents at the Fermi surface at each 
Q even allowing for Fermi surface fluctuations, which may be integrably singular at low 
frequencies. 

However, this does not imply that, in the presence of interactions, d a and p a (or Afcpg-) 
remain the appropriate canonically conjugate variables. These are variables which measure, 
respectively, the particle number at a particular patch on the Fermi surface and given spin 
and the phase of the wave function at the Fermi surface. If there is a finite phase shift 
for forward scattering of opposite spin electrons, as we have shown [jl], |7j, the order of doing 
these operations matters. If we add a particle of up spin, the phase of the down spin wave 
function depends on whether the particle of up spin was added before or after the phase 
was measured. The failure of commutation for opposite spins is the phase shift rj/ir, just as 
adding a particle of up spin below the Fermi surface enforces a change in up-spin phase by 
the amount 7r. We may express this by writing the free particle commutator in matrix form: 



if)(x) = pe 



,W{x) 



(4) 
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while 

[Pa, d a ] interacting = «7T I ^ " J X 6(r - r') S(Q - Q'). (6) 

Let us explain these equations in detail. (5) means in the one-dimensional model that if 
we insert an extra particle into the Fermi sea at a point r, because of the exclusion principle 
the wave-function at the Fermi surface (which is the basic interpretation of equation (4)) 
must have an extra node inserted into it near r, hence the phase difference between left and 
right-going (or in — and out — going) waves must shift by n as a consequence. Hence after 
we insert one particle in p, 9 will change by n, but not vice versa: one is the generator of 
displacements of the other. (6) must be interpreted in exactly the same way. The insertion 
of an up-spin particle at r, near Q, means that the phase of the down spin wave at Q is 
shifted by rj, while the up-spin wave is shifted by ir. This means that pj and 0±, p\ are no 
longer canonically conjugate; the correct canonically conjugate variables are proportional to 

and 

The equations of motion of the charge and spin bosons follow from the commutation 
relations and the Hamiltonian, which as we explained is simply the original kinetic energy, 
the interaction terms being completely subsumed in the C.R. The Hamiltonian is as for free 
particles the one given by Haldane, 

n = ^Jdnjd D rv F (n)[Ak F (n),r,t)] 2 
= UdnJ2v F ip 2 (q,n) + q 2 9 2 (q,n)]. (9) 

q 

Then 

[H,e c , a {q,n)\=v c , aq e c , a {q,n) (10) 

with 

«. = Mn)(i-— ) 
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v e = v F (n)(i + ^) (n) 

f qv c (fl) 

and bosons are left as harmonic oscillator variables with frequencies < 

{ qv s (fi) 

For free particles the Fermion operator is made up from bosons via the formula 

Free: ^(r) = Po e-^ {dc+<jds) (12) 
which gives the Green's function 

Gfree = j 1 = j 1 = (t>c = V a ) (13) 

^Jr -v s t + i/A ^r-v c t + ^ 

But we cannot assume that the connection between interacting electrons and the modified 
bosons obeys (12). The coupling of the two Fermi surfaces which leads to the modified 
C.R. means that (12) creates an object which can be thought of as a "pseudoelectron" with 
the suitable backflow caused by the fractional opposite-spin hole which accompanies it, so 
it describes an exact eigen-excitation of electron-like character moving in the exact ground 
state. These excitations are analogous to bosonized versions of the exact eigen-excitations of 
charge (/j) and spin ( J a ) of the Lieb-Wu solution of the ID Hubbard model. (The discussion 
here was foreshadowed in Y. Ren's thesis |T2]| .) These ladders of excitations can be described 
in terms of appropriate bosons since they have linear energy-momentum relations near zero 
energy, and these are the bosons which we have derived. But the actual electron operator 
creates a physical electron, not the pseudo-electrons described by these bosons, and hence 
must have the backflow compensated out. This leads to the fractional exponents in the 
Green's function and other correlation functions characteristic of the Luttinger liquid. As 
in the ID case (as shown in Ren's thesis) the coefficients may be deduced from conservation 
laws and from the Luttinger theorem of incompressibility of the Fermi sea in momentum 
space. 

Note that the "pseudo-electron" has the quantum numbers of a true electron, and in 
fact it is one of the packet of exact eigenstates created when a true electron is inserted 
at the appropriate momentum, though with vanishing amplitude as L — > oo. When a real 
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electron is added, a cloud of particle-hole excitations in addition to the two semions is excited, 
analogous to the cloud of particle-hole excitations which causes the x-ray edge anomaly. This 
is the "backflow". The modified commutation relations of the charge and spin bosons still 
leave them as a bosonic description of particles which are "semions" in the sense that two 
of them make an electron. The transformation which diagonalizes the C.R. is not modified 
from the free particle case, i.e., it is independent of rj. 

This is essentially because we maintain Luttinger's theorem of incompressibility as a 
constraint, so that no net down-spin particles are removed by the scattering process: they 
are merely redistributed in momentum space, which is the "backflow" we must now calculate. 
^ particles are displaced from the neighbourhood of the scatterer particle at k f, and we 
must find how they displace the Fermi surface bosons, i.e., how the phases are shifted at the 
Fermi surface. But first we must take into account some consequences of the non-Abelian 
spin symmetry which we have been ignoring so far. 

A key theorem of the bosonization technique follows from the symmetry properties of 
the states at the Fermi surface. As we said before, the existence of a Fermi surface implies 
separate conservation of each component of spin at each point on the Fermi surface. But 
spin conservation must remain independently of the choice of axes, and we must be able 
to choose the axes at each point independently. A related requirement is that Kramers 
degeneracy of the spin at each point of the Fermi surface independently must be maintained. 
This is not possible if spin at different Fermi points is coupled relevantly as uj — > 0. As is seen 
in the ID Hubbard model, this implies that the spin bosons cannot acquire an anomalous 
dimension, and must retain the same semionic character that they have for free Fermions. 
In our situation, this expresses itself by the observation that our scattering calculation is 
slightly incomplete. We have not required formal spin rotation invariance (SU(2) symmetry) 
of the S-matrix for scattering, which requires that the phase-shift have the form 



and allows for a spin-flip scattering, which we have so far ignored, of half the magnitude r] of 
the potential term. This requires the scattering to take place entirely in the singlet channel, 




(14) 
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rather than the up-down channel, as we have implied in our discussion so far. Our previous 
picture left us with one spin fcj plus a hole of magnitude 2 in fcj. This left 1 + a f spins, but 
now we have spin-flip scattering of 77/2 giving ^- missing down spins and 1 — ^ up spins 
or one net spin. Correspondingly, this gives matching currents of up spin in the scattered 
channels which leaves us with displacements only of charge, not spin, bosons in the backflow. 
The co-moving hole of magnitude rj/ir is now in the charge channel. 

In the actual ID Hubbard model, this theorem is satisfied only to logarithmic accuracy, 
leading to (lntu) -1 and (lng) -1 corrections to power laws; the relevant coupling constant 
goes to zero only logarithmically. We expect the same pathology in 2D. But dominant 
power laws will be correctly determined by bosonization. (All of this was foreshadowed in 
Haldane's "Luttinger Liquid" treatment of the ID Hubbard model. [0) When the spin-flip 
component is taken into account, we now can determine how the phases at the Fermi surface 
are shifted, specifically when we insert an electron at Q, q in order to calculate the one- 
particle Green's function. The rule is very simple: we calculate the phase shifts we would 
have expected using naive up-spin down-spin scattering, and replace these by phase-shifts 
in the pure charge channel. Let us first discuss the ID case, which was worked out by Ren 
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In ID, the amount of charge ^ which is displaced from the state k = kp — q appears, 
half at the left-hand Fermi point and half at the right, i.e., 7?- at each. These components 
multiply the Green's function by the factor 



which gives, in space-time representation, a factor 

1 



-(x 2 -v 2 c t 2 )i ^ 

which, has the maximum exponent (|) 2 x | = ^, as pointed out by Ren. This gives the 
famous Fermi surface smearing exponent 2 x 1/16 = 1/8, in the strong coupling case and 
with the strictly local interaction appropriate to the Hubbard model. 
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These two displacements are the total backflow. The net momentum of the backflow is 
zero, and the net charge i]/n , as it must be. 

The situation in 2D is not quite so simple. Again, we recognize that rj/ir worth of 
charge boson — i.e., i]/7r enclosed by an "internal Fermi surface". — has been displaced from 
the region of momentum k. We may calculate the displacement of a circular Fermi surface 
which would result from elastic incompressible deformation of the lattice of /c-values. (We 
use a circular FS for illustrative purposes.) This would give us 

s%> = k ' lk . A_ 

{k'-kf 2tt 2 

and 

5kf) M^^^ (16) 

See Fig. (1). If k is chosen at 9 — 0, and k = k F — e 



., ,, v ek F ri ii 1 — cos# i) r , M 11 



In this case, half of the displacement is in the forward direction, and half is a uniform displace- 
ment of the Fermi level — essentially an s-wave, equivalent to isotropic potential scattering. 
This, however, is not quite the whole story. In one dimension the backflow compensated 
the charge and the momentum exactly, since the left-and right-moving pieces were identical. 
Here, however, we have an uncompensated momentum of the forward-moving wave, j- x kp. 
The correct displacement satisfying the Luttinger-Ward theorems is not merely a dilation of 
the momentum lattice, but a rigid displacement of — ^ kp as well. 

The simple uncompressible dilation of the Fermi surface which we postulated in (16) is 
too simple: the interactions must satisfy momentum as well as particle conservation and 
so the backflow must carry no net momentum, as in ID. The relative s-wave channel must 
carry momentum — j^kp which compensates the extra momentum of the 5-function peak 
at kp. This is equivalent to a uniform translation of the Fermi surfaces, which is a simple 
unitary transformation (multiplication of all states by a common factor) and does not lead 
to any anomalous dimensions. On the other hand, the s-wave dilation does do so, and the 

10 



anomalous dimension of the Green's function is, as in ID, {^) 2 x | = a, < a < 1/8. 
Another way of describing this part of the backfiow is as a Fermi surface shift proportional 
to (1 — 2 cos#) rather than simply to 1. This is not a scattering in the p-wave channel, rather 
it is more like a "Mossbauer" zero-phonon, coherent recoil of the Fermi whole. 
The form of the Green's function is quite different from ID: it will look something like: 

n(Q) is the Fermi surface normal unit vector at Q, and cosd = h(Q)r. Stationary phase will 
ensure that G(r,t) comes almost entirely from the "patch" n(fi)||r. 

Experimentally, several hints suggest that a > 1/8 in fact, in the cuprates. We must 
not be surprised by the parallel-spin interaction also being finite and repulsive, which will 
enhance the charge-channel backfiow without affecting spin properties except to lower v s 
further, and make the electrons even less Fermi-liquid. For the Hubbard model there is a 
fixed relation between r] c and i] s in (14), but in the physical case rj can be larger. 

Most of the physical phenomena which depend on G and other correlation functions can 
be calculated using the simple homogeneity property: 

G=4r- H~) • (2°) 

This determines the infrared spectrum in parallel and perpendicular polarizations |TJ], [TJJ, 
and the Fermi surface smearing; a similar property will give the exponent for \fT\. Only 
ARPES requires the full G, and this will depend critically on details of the single-particle 
dispersion and Fermi surface, so will require a separate investigation. 

With (19) we have in principle the asymptotic solution of the 2D electron gas with a local, 
repulsive interaction. This is expected to be valid in the regions of the phase diagram of the 
Hubbard model reasonably far from half-filling (where umklapp terms are important and 
can pin down the charge bosons) and U — > oo at high density, where ferromagnetic coupling 
of Landau mean field type will possibly be important, and lead to Nagaoka ferromagnetism. 
Finally, we exclude strong magnetic fields, "strong" being enough to allow interference after 
a full cyclotron orbit, i.e., we require u c r << 1. Under this condition transverse gauge 
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transformations are simple reparametrizations of the Fermi surface and meaningless; i.e., 
the Fermi surface and anyons are mutually incompatible. fl c r >> 1 destroys the symmetries 
implicit in the Fermi surface and causes gaps in the spectrum which are incompatible with 
bosonization. With Q c r < 1 bosonization is the gauge theory of the interacting Fermi system; 
there is no meaningful other. 

Khveshchenko has argued that in > 2 dimensions the equations of motion of the bosons 
are a very crude approximation valid only for very small q and uj. This is clearly so in our 
approach, since the "5-function" in equation (17) is actually of width q. We have argued that 
"Chern-Simons" types of terms are not important if uj c t < 1, but in so far as charge and spin 
velocities differ these can be effects such as those we have postulated in the past caused by 
mixing of bosons over a finite area of the Fermi surface, when electrons of finite q = k — Icf 
are excited. Thus the above is a first approximation to a much more complex theory which 
we do not yet have under control. Nonetheless it seems the only way to proceed. 
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